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Abstract 

Events stemming from the pair-production of Z bosons in e"'"e~ collisions are 
studied using 217.4 pb^^ of data collected with the L3 detector at centre-of-mass 
energies from 200 GeV up to 209 GeV. The special case of events with b quarks is 
also investigated. 

Combining these events with those collected at lower centre-of-mass energies, the 
Standard Model predictions for the production mechanism are verified. In addition, 
limits are set on anomalous couplings of neutral gauge bosons and on effects of extra 
space dimensions. 
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1 Introduction 



The pair-production of Z bosons in e"'"e~ collisions at LEP was observed [1] with the L3 de- 
tector [2] once the accelerator centre-of-mass energy, ^/s, exceeded the production threshold 
of 2mz, where mi denotes the Z boson mass. Numerous studies from data samples collected 
at the steadily increasing ^/s and integrated luminosities were reported by L3 [3, 4] and other 

collaborations [5]. 

In the Standard Model of the electroweak interactions [6], the Z pair-production is described 
at the lowest order by two i-channel Feynman diagrams with an internal electron^^ leg, collec- 
tively denoted as NC02. A wider definition is used in this Letter: all diagrams leading to two 
fermion-antifermion pairs are considered and kinematic restrictions which enhance the NC02 
contribution are enforced. Results in the NC02 framework are also given. 

The study of Z pair-production offers a further test of the Standard Model in the neutral 
gauge boson sector and is of particular relevance as this process constitutes an irreducible 
background in the search for the Standard Model Higgs boson at LEP. Events with Z boson 
decaying into b quarks have a signature similar to those originated by the process e"'"e~ — > ZH — > 
ffbb. Their selection and the measurement of their cross section validate the experimental 
procedures used in the search of the Standard Model Higgs boson. 

Z pair-production allows the investigation of the anomalous triple neutral gauge boson 
couphngs ZZZ and ZZ7 [7] , forbidden at tree level in the Standard Model and tests new theories 
like possible effects of extra space dimensions [8]. 

This Letter describes measurements at two average values of ^/s, 204.8 GeV and 206.6 GeV, 
corresponding to integrated luminosities of 78.5 pb~^ and 138.9 pb~^, respectively. Hereafter, 
these data samples are denoted as the 205 GeV and 207 GeV energy bins. Gombined results 
from the full Z pair-production sample collected with the L3 detector are also given in the 
comparison with Standard Model expectations and for constraints on New Physics. 

2 Monte Carlo simulations 

The EXCALIBUR [9] Monte Garlo program is used to model signal and background neutral- 
current four-fermion processes. The Z pair-production process is defined as the subset of the 
four-fermion generated phase space satisfying the following kinematics cuts [1,3,4]. The in- 
variant mass of fermion-antifermion pairs is required to be between 70 GeV and 105 GeV. For 
events with two identical pairs, at least one of the possible pairings has to satisfy this condition. 
For the uddfi, esse and Vii'^Vfi^ [i = e, /x, r) final states, the masses of the pairs that could 
originate from a W decay have to be either below 75 GeV or above 85 GeV. The polar angle 
9e of generated electrons is required to satisfy | cos6'e| < 0.95. 

The Z pair-production cross section is calculated to be 1.07 pb and 1.08 pb for the 205 GeV 
and 207 GeV energy bins, respectively. Following a comparison with the GRC4F [10] Monte 
Carlo generator, and taking into account the modelling of initial state radiation, an uncertainty 
of 2% is assigned to these calculations. 

The cross section for final states with b quark pairs is significantly smaller. Gombining the 
two energy bins a cross-section of 0.30 pb, also with an uncertainty of 2%, is expected for an 
average centre-of-mass energy ^/s — 205.9 GeV. 

^^In this Letter, the word electron is used for both electrons and positrons. 
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Four-fermion events generated with EXCALIBUR which do not satisfy the signal defini- 
tion are considered as background. Background from fcrmion pair-production is described by 
KK2f [11] for the processes e'^e^ qq(7)i c+e~ fi'^/j^{j) and e+e" — > t~^t~ {'-/), and BH- 
WIDE [12] for e"'"e~ — > e"'"e~(7). Background from charged-current four-fermion processes is 
generated with EXCALIBUR for the ei/gqq' final state and with KORALW [13] for W pair- 
production and decay in final states not covered by the simulations listed above. Hadron and 
lepton production in two-photon processes is modelled by PHOJET [14] and DIAG36 [15], 
respectively. 

The L3 detector response is simulated using the GEANT program [16], which takes into ac- 
count the effects of energy loss, multiple scattering and showering in the detector. GHEISHA [17] 
is used for the simulation of hadronic interactions. Time dependent detector inefficiencies, as 
monitored during the data taking period, are also reproduced. 

3 Event selection 

All visible final states of Z pair-production are investigated. For the qquu, i'^£~i'u and 
t^l~£''^£'~ final states, criteria are used which are similar to those developed at = 189 GeV [3] 
and = 192 — 202 GeV [4]. For the qqq'q' and qq£+£~ final states, improved analyses are 
devised. All selections rely on the identification of two fermion pairs with masses close to mz- 

Electrons are identified by requiring a well isolated electromagnetic cluster in the electro- 
magnetic calorimeter with an associated track in the tracking chamber. To increase efficiency, 
the track matching requirement is relaxed in some selections. 

Muons are reconstructed from the coincidence of tracks in the muon spectrometer and the 
central tracker which are in time with the beam crossing. Energy depositions in the calorimeters 
which are compatible with a minimum ionising particle (MIP) and have an associated track in 
the central tracker are also accepted. 

Taus are identified by their decays either into electrons or muons, or into hadrons detected 
as narrow and isolated low multiplicity jets associated with one, two or three tracks. 

Quark fragmentation and hadronisation yields a high multiplicity of calorimetric clusters 
and charged tracks. These are grouped into jets by means of the DURHAM algorithm [18]. 
The number of reconstructed jets depends on the thresholds ymn for which a m-jet event is 
reconstructed as a n-jet one. 

The tagging of b quarks [19] relies on the reconstruction of the decay vertices of weakly 
decaying b-hadrons with the silicon vertex detector and the central tracker. The shape and 
particle content of the associated jets are also considered. 

The four-momenta of neutrinos are derived from all other particles measured in the event, 
making use of the hermeticity of the detector. 

3.1 The qqq'q' channel 

The study of the qqq'q' channel [20] starts by selecting high multiplicity events with a visible 
energy, E^ig^ satisfying 0.75 < E^isl \fs < 1.35. The energy imbalance in the directions parallel 
and perpendicular to the beam axis have to be less than Q.2Ems and 0.25£'^is, respectively. 
These criteria suppress fermion-pair production, two-photon interactions and four-fermion final 
states with leptons. To further reduce boson pair-production with leptons, events with isolated 
electrons or muons with an energy larger than 40 GeV are rejected. Events with isolated 
photons of energies above 25 GeV are also rejected. 
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The remaining events are forced into four jets. A kinematic fit which imposes four- momentum 
conservation is performed to improve the di-jet mass resolution. Among the three possible jet 
pairings, the pairing i is chosen which minimizes: 

where Sj and Aj are the di-jet mass sum and differences and o"!^^ and cr]^^^ their resolutions, 
determined from Monte Carlo. Only events for which > 165 GeV are retained. 

The remaining background is formed by events from the e+e" — > W+W" and e+e~ — > qq(7) 
processes. A likelihood, I>fz , is built which combines ten variables: the event sphericity, the 
largest triple jet boost [21], the largest jet energy and boost, the largest energy difference 
between any two jets, the opening angle between the most and least energetic jets, log 7/34, the 
mass M5C from a kinematic fit with equal mass constraint, the absolute value of the cosine 
of the polar angle of the event thrust vector, | cos 9t\, Xzz corresponding variable for 

the W-pair hypothesis, Xww- The distributions of Lfl' for data and Monte Carlo are shown 
in Figure 1. Events with Lfl' < 0.1 are mostly due to the e+e" qq(7) process and are not 
considered in the following. 

A second likelihood, Lzz, is built to further exploit the difference between Z pair-production 
and the residual background from W pair- production. It uses seven variables: Lfl' , Ezz, the 
corresponding variable for the W-pair hypothesis, Eww, the three to four-jet threshold for the 
JADE clustering algorithm [22], the event thrust and the value of the b-tag variable for the 
two jets with the highest probability to originate from b quarks. 

The distributions of Lzz for data and Monte Carlo are shown in Figure 2a. Table 1 lists 
the yield of this selection for Lzz > 0.2, which corresponds to an efficiency of 55.4%. 

3.2 The qquP channel 

The selection of the qqui? channel is identical to that performed at ^/s — 192 — 202 GeV [4]. 
High multiplicity events with large missing energy and momentum and no high energy isolated 
electrons, muons or photons are selected. They are forced into two jets and a constrained fit 
which enforces the hypothesis that the missing four-momentum is due to a Z boson is applied. 

Finally, an artificial neural network singles out Z pair-production events from background. 
Its input variables include event shape variables that differentiate a two-jet from a three- 
jet topology, the sum of visible and missing masses, the masses of the two jets, the missing 
momentum and the energy in a 25° azimuthal sector around the missing momentum vector. 

Figure 2b shows the output of the neural network, NNout, for data and Monte Carlo. The 
results of this selection are summarised in Table 1 for a benchmark cut NNout > 0.5, which 
corresponds to an efficiency of 46.2%. 

3.3 The qq£+£- channel 

The study of the qqe'^e", qqAi"*"/^" and qqr+r" final states [23] proceeds from a sample of high 
multiphcity events, well balanced in the planes parallel and transverse to the beam direction. 
Background from two- photon interactions is rejected by requiring | cos^tI < 0.98. Events from 

the e+e" qq(7) process with hard initial state radiation photons are reduced by requiring 
the effective centre-of-mass energy [24], V^, to be greater than 0.55y^. Remaining two-jet 
events are suppressed by a cut on the event thrust. 
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To fully reconstruct the qq£'^i~ final state, after identifying an electron, muon or tau pair 
in the event, the remaining clusters are forced into two jets. A kinematic fit which enforces 
energy and momentum conservation and equal mass, M5C, of the hadronic and leptonic systems 
is performed. To cope with different background contributions, different selection criteria are 
applied for the three final states. 

The qqe"'"e~ selection requires low transverse energy imbalance and a sum of the energies of 
the two electrons close to The variable (£^1 + E2 — E3 — E4)/{Ei + E2 + E3 + E4), is also 

considered, where Ei denotes the decreasingly ordered jet and lepton energies. This variable has 
low values for the signal, where the energy is uniformly distributed among the four particles, 
and large values for the background from the e+e~ —>■ qq{j) process. An efficiency of 73.9% is 
reached. 

The qq/i+//~ final state has a low background contamination, almost entirely rejected by 
requiring a large energy for the lowest energetic muon, expected to be soft for background 
events. This selection has an efficiency of 60.4%. 

The qqr'^T" selection is affected by a larger background. It requires the invariant mass of 
the hadronic system prior to the kinematic fit to be compatible with mz and a large rest frame 
angle between the taus. In addition, the sum of the di-jet and di-tau masses after the kinematic 
fit has to be close to 2mz. This selection accepts 28.0% of the qqr+r" final states, as well as 
2.4% and 4.8% of the qqe+e^ and qqn^n~ final states, respectively. 

Table 1 presents the combined yield of all selections, whose overall efficiency is 55.4%. 
Figure 2c shows the M^c distributions for data and Monte Carlo. 

3.4 The £+£.-vD channel 

Only final states with electrons and muons are considered with a selection identical to that 
performed at lower ^/s [4] . A pair of acoplanar leptons is selected in low multiplicity events with 
large missing momentum pointing away from the beam axis. To improve efficiency, electrons 
are not required to have an associated track. To reduce the background, no MIP candidates 
are accepted in the muon selection. Both the lepton visible mass, M^, and recoil mass, Mrec, 
must be consistent with mz- These criteria suppress background from fermion pair- production. 
Residual background from four- fermion processes is reduced by performing a fit which constrains 
the leptons to originate from a Z boson and requiring the recoil mass to be close to mz. 

For signal events, the sum M« + Mrec should be close to 2mz. The distributions of this 
variable for data and Monte Carlo, combined with results from the £^£^£'^£'^ channel, are 
shown in Figure 2d. An efficiency of 25.3% is achieved and the results of the selection are 
reported in Table 1. 

3.5 The £+£-£'+£'- channel 

The selection for the £'^£~£!^£!~ channel aims to retain a high efficiency to compensate for the 
low branching ratio. The same criteria of Reference 4 are applied to select low multiplicity 
events with four or more loosely identified leptons, with energy above 3 GeV. Events must 
contain at least one electron or muon pair. Electrons without an associated track are accepted 
while MIPs are not considered to form these pairs. 

The lepton pair with mass closest to mz is selected and both M(i and M^ec are required 
to be compatible with mz. Background from fermion pair-production is reduced by imposing 
upper bounds on the opening angle of the leptons of this pair. 
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The sum + M^ec is used as a final discriminating variable. Its data and Monte Carlo 
distributions are presented in Figure 2d, together with those from the i^i^uu channel. The 
yield of the selection is summarized in Table 1 and corresponds to an efficiency of 33.6%. 

4 Results 

4.1 Measurement of the ZZ cross section 

The cross sections for each energy bin and each final state are derived [4] with a fit to the 
final discriminating variables and are presented in Table 1 together with the Standard Model 
predictions. In the presence of fluctuations for channels with low statistics, an upper limit [4] 
on the cross section is given. Fixing the relative contributions of all channels to the Standard 
Model expectations, the Z pair-production cross section is extracted and also presented in 
Table 1. All the measured cross-sections agree with their Standard Model predictions. 

4.2 Study of systematic uncertainties 

Several sources of systematic uncertainty are considered [4] and listed in Table 2. Systematic 
effects correlated among channels arise from uncertainties on the detector energy scales, on the 
signal modelling, as derived from a comparison between EXCALIBUR and GRC4F and on the 
prediction of the background level. This is studied by varying the expected cross sections for W 
pair-production, jet production, the eu^qq' and four-fermion processes by 0.5%, 5%, 10% and 
5%, respectively. The qqq'q' channel is affected by uncertainties on the charge multiplicity and 
the simulation of the b-tag discriminant. Other sources of systematic uncertainty, uncorrelated 
among the channels, are the signal and background Monte Carlo statistics, detailed in Table 3 
and the accuracy of the simulation of the selection variables and of those used for the lepton 
identiflcation. 

Including all systematic uncertainties, the measured cross sections read:^^ 



Figure 3 presents these values together with lower energy measurements [1,3,4]. 
4.3 Final states with b quarks 

The final discriminant of the qqq'q' analysis, plotted in Figure 2a, shows a high sensitivity 
to final states containing b quarks. In order to tag hhuu and hh£^fr final states, the b-tag 
information of the two jets are combined [4] to form the discriminating variables shown in 
Figure 4. For the hhuu final state, the value of the variable NNgut is also considered in the 
combination. The cross sections for Z pair-production with b quarks in the final states are 

^^In the NC02 framework, the cross sections are derived as: 



where the Standard Model expectations, consistent among the EXCALIBUR, ZZTO [25] and YFSZZ [26] 
programs, are assigned an uncertainty of 2%. 



(Tzz(205 GeV) 
(7zz(207GeV) 



0.78 ± 0.20 (stat.) ± 0.05 (syst.) (SM : 1.07 ± 0.02 pb) 
1.10 ± 0.17 (stat.) ±0.07 (syst.) (SM : 1.08 ± 0.02 pb). 



(7^^°2(205GeV) = 0.77 ±0.20 (stat.) ±0.05 (syst.) (SM : 1.05 ± 0.02 pb) 
(7^1^02 (207 GeV) = 1.09 ± 0.17 (stat.) ± 0.07 (syst.) (SM : 1.07 ± 0.02 pb) 
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determined from a fit to these variables and listed in Table 4. Their combination gives a total 
cross section: 

t^zz^bbx (205 - 207 GeV) = 0.24 ± 0.09 (stat) ± 0.03 (sys). 
The systematic uncertainty follows from the sources discussed above and is detailed in Table 2. 

4.4 Combined results 

The ratio between measured and expected cross sections, i?zz = a^'^^ ja^^ ^ is calculated including 
lower energy data [1,3,4] as: 

i?zz (183 - 209 GeV) = 0.93 ± 0.08 (stat) ± 0.06 (sys) . 

Systematic uncertainties include correlations among different data samples. The predictions 
are in agreement with the measurements with a precision of 11%. 

Figure 5a shows the distribution of the reconstructed mass Mi of the Z boson, and Figure 5b 
the absolute value of the cosine of the observed production angle for the full Z pair-production 
sample. The cuts L^i > 0.85 and NNout > 0.8 are applied to data decsribed in this Letter. 
Data at lower energies [1,3,4] are also included. 

5 Limits on physics beyond the Standard Model 

5.1 Anomalous couplings 

Assuming on-shell production of the two Z bosons, anomalous ZZV couplings are parametrised [7] 
by the coefficients f^, with i = 4, 5 and V = 7, Z. The // coefficients correspond to CP viola- 
tion and the ones to CP conservation. All fy coefficients are zero in the Standard Model. 
Each event of the signal Monte Carlo distributions presented in Figure 2 is reweighted [3] to 
simulate anomalous values of the coefficients. The full phase space of the Z boson pair, 
as reconstructed from the jet and lepton four-momenta, is used. A fit to these distribution is 
performed leaving one coupling free at a time and fixing the others to zero, yielding the 95% 
confidence level limits: 

-0.48 < /I < 0.46; -0.36 < < 1.03; -0.28 < f2 < 0.28; -0.40 < < 0.47, 

compatible with the Standard Model expectations. Lower energy data [1,3,4] are included in 
the fit. Figure 6 presents results of simultaneous fits to couphngs with the same CP eigenvalue. 

5.2 Extra space dimensions 

A recent theory [27] , dubbed "Low Scale Gravity" , proposes a solution to the hierarchy prob- 
lem by postulating a scale Ms for the gravitational interactions which is of the order of the 
electroweak scale. Extra space dimensions are a consequence of this theory. In this scenario, 
spin-two gravitons contribute to the Z pair-production [8], interfering with the Standard Model 
production mechanism. The Z pair-production cross sections presented in this Letter and those 
measured at lower energies [1,3,4] are fit with a combination of Low Scale Gravity and Standard 
Model contributions. A lower 95% confidence level limit on the scale M5 of 0.7 TeV is obtained. 
It holds for both constructive and destructive interference between the Low Scale Gravity and 
the Standard Model contributions. 
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Channel 


x/s CGeV) 








af^^ fpb) 


a^^ fpb) 


noo'o' 

M.M.M. M. 


205 


166 


24.9 ± 0.0 


140.5 ± 0.4 


38+^-?? 


0.51 




207 


300 


46.6 ± 0.0 


255.1 ± 0.6 


0.35 r\ 1 /I 
— U.i4 


0.52 




205 


13 


10.8 ± 0.1 


11.0 ± 0.1 


< 0.24 


0.30 




207 


36 


19.3 ± 0.1 


18.7 ± 0.1 




0.30 




205 


10 


6.8 ± 0.1 


1.6 ± 0.1 




0.16 




207 


18 


12.3 ± 0.0 


3.0 ± 0.1 


0.19+n^ 


0.16 




205 


2 


0.9 ± 0.0 


0.6 ± 0.0 


o.o8l°:°i 


0.04 




207 


3 


1.6 ± 0.1 


1.2 ± 0.0 


0.05lH^ 


0.04 


£+£-£'+£'- 


205 





0.4 ± 0.0 


0.2 ± 0.0 


< 0.11 


0.02 




207 


1 


0.7 ± 0.0 


0.4 ± 0.0 


o.o8l°:°i 


0.02 


e+e- ^ ZZ 


205 


191 


43.8 ± 0.1 


154.0 ± 0.4 


0.78 ±0.20 


1.07 




207 


358 


80.4 ± 0.1 


278.4 ± 0.6 


1.10 ±0.17 


1.08 



Table 1: Number of observed data events, A^/j and signal, Ns, and background, Nb-, expected 
Monte Carlo events in the two energy bins. The benchmark criteria L^z > 0.2 and NNout > 0.5 
are applied for the qqq'q' and qquu final states, respectively. Uncertainties are due to Monte 
Carlo statistics. Measured, cr^'*, and expected, a*^, cross sections are also given. Limits on cr-^'* 
are at the 95% confidence level. 





Sazz (%) 


(^CTzZ-^bbX (%) 


Correlated sources 






Energy scale 


3.1 


2.4 


Theory predictions 


2.0 


2.0 


WW cross section 


0.4 


0.5 


Four-jet rate 


1.4 


2.7 


Wez/ cross section 


1.1 


0.8 


Four-fermion cross section 


0.5 


0.5 


Uncorrelated sources 






Charge multiplicity 


1.3 


2.3 


B-tag 


2.5 


11.3 


Monte Carlo statistics 


1.9 


3.1 


Simulation 


3.5 


3.5 


Total 


6.4 


13.2 



Table 2: Systematic uncertainties on azz and criz^hhx- 
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qqi+i- 




qqq'q' 




i+i-i'+i'- 


Signal MC statistics (ctzz) 


0.5% 


0.4% 


<0.1% 


3.0% 


1.0% 


Background MC statistics (o"zz) 


3.4% 


0.5% 


0.2% 


0.1% 


0.2% 


Signal MC statistics (o"zz^bbx) 


3.9% 


1.8% 


2.1% 






Background MC statistics (o"zz^bbx) 


3.4% 


1.4% 


2.2% 






Simulation 


4.8% 


3.2% 


1.6% 


0.8% 


1.8% 



Table 3: Sources of uncorrelated systematic uncertainties on ctzz and cTzz-^hhx- 







hhuu 


qqbb 


Measured cross section (pb) 
Expected cross section (pb) 


0.032 ± 0.027 
0.035 


< 0.108 
0.065 


0.185 ± 0.074 
0.201 



Table 4: Cross sections for final states with b quarks. The limit is at 95% confidence level. 
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Figure 1: Distribution of the likelihood Lfl' used for the qqq'q' selection. The signal and 
background Monte Carlo distributions are normalised to the expected cross sections. 
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Figure 2: Distribution of the final variables used for the measurement of the cross-section for 
the a) qqq'q' b) qqz/i/ and c) qq£"'"£~ final states. The sum of the t^i~vv and final 
states is given in d). The signal and background Monte Carlo distributions are normalised to 
the expected cross sections. 
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Figure 3: Measurements and predictions for the e"'"e —>■ TIL and e"^e TIL - 
sections as a function of ^/s. An uncertainty of 2% is assigned to the predictions. 
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Figure 5: Distributions in data and Monte Carlo at all LEP centre-of-mass energies above 
the Z pair production threshold of a) the reconstructed mass M and b) the absolute value of 
the cosine of the production angle dz- Cuts on the qqq'q' final discriminant and on the qquu 
neural network output are applied. The signal and background Monte Carlo distributions are 
normalised to the expected cross sections. 
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Figure 6: Results of a simultaneous fit to anomalous coupling parameters with the same CP 
eigenvalue. The Standard Model (SM) expectations are also indicated. 
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